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ABSTRACT 

A large number of massive stars are known to rotate rapidly, resulting in significant 
distortion and variation in surface temperature from the pole to the equator. Radia- 
tively driven mass loss is temperature dependent, so rapid rotation produces variation 
in mass loss and angular momentum loss rates across the surface of the star, which is 
expected to affect the evolution of rapidly rotating massive stars. In this work, we use 
ZAMS stellar models to investigate the two dimensional effects of rotation on stellar 
mass loss, using two common prescriptions for radiatively driven mass loss. The associ- 
ated loss of angular momentum from these models is also considered. Using 2D stellar 
models, which give the variation in surface parameters as a function of co-latitude, 
we implement two different mass loss prescriptions describing radiatively driven mass 
loss. We find a significant variation in mass loss rates and angular momentum loss as 
a function of co-latitude. We find that the mass loss rate decreases as the rotation 
rate increases for models at constant initial mass, and derive scaling relations based 
on these models. When comparing 2D to ID mass loss rates, we find that although the 
total angle integrated mass loss does not differ significantly, the 2D models predict less 
mass loss from the equator and more mass loss fr om the pole than the ID predictions 
using von Zeipel's law. As a result, rotating models lose less angular momentum in 
2D than in ID, which will change the subsequent evolution of the star. The evolution 
of these models will be investigated in future work. 
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1 INTRODUCTION 

Radiatively driven mass loss plays an important role in de- 
termining the properties of the circumstellar environment, 
and is also expected to have a significant impact on the sub- 
sequent evolution of the star itself. Predominantly impor- 
tant in massive stars, the radiatively driven stellar winds 
are an important factor in determining the course of evo- 
lution and final end state of the star. However, many mas- 
sive stars are known to be rapidly rotating, which produces 
variation in the flux and effective temperature over the stel- 
lar surface. This variation in turn, will affect the tempera- 
ture dependent mass loss rates, producing lati tudinally de- 
pendent stellar winds (see, for example lOwocki et al.iil998: 
iDwarkadas fc Owockill2002l '). The effects of rotation on stel- 
lar mass loss rates can greatly influence t he later evolution of 
a sta r through angular momentum loss l|Mevnet fc Maedeil 
L2OOO'). It is therefore important for us to understand how ex- 
actly mass loss rates depend on rotation, and the subsequent 
effects on evolution. 

In general, the effects of stellar rotation on wind pro- 
files are not well understood. Research into this prob- 
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lem is still in its infancy, and the results can be 
highly contested. Several studies have shown that at very 
high rotation rates, the wind can switch to a slower 
denser equ a torial outflow JPelupessv. Lamers fc Vinkll2000l: 



Curd I2OO4I: ICure fc Riall l2004l: ICure. Rial, fc Cidald 



20051 : 



Madura. Owocki fc Feldmeieijl2007l ). which may explain the 
observed disks around B[e] stars. However, other groups find 
that when the effects are modelled in 2D, the opposite is 
true, an d the wind density is higher at the pole than the 
equator (|Petrenz fc PulsllioOOl '). However, in this work, we 
focus instead on how rotation affects the mass loss rates 
from the surface of the star, and the implications these mass 
loss rates have for the structure and evolution of the star, 
rather than the winds. 

Early investigations into latitudinally dependent mass 
loss suggested that the reduction in effective gravity near the 
equator would produce an enh ancement in the equatorial 
density l|Friend fc Abbofr3ll986l ). The increase in mass lost 
at the equator would carry away more angular momentum, 
and the star would spin down faster. However, these results 
assumed the flux to be constant across the surface of the 
star. When the effects of gravity darkening are taken into 
account, (using, e.g., von Zeipel's law,) the opposite was 
found to be true, and the star was found to lose most of it's 
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mass from the polar regions l|Dwarkadas fc Owocki|[2"002l ) . In 
this case, the star will spin down more slowly, which could 
lead to relatively rapidly rotating stars in later phases of 
evolution, such as the luminous blue variable (LBV) or Wolf- 
Rayet (WR) phases. 

A change in the mass loss rate will change the amount 
of angular momentum lost from a star, causing it to spi n 
down differently. As discussed bv lMevnet fc Maederl (l2000l l. 
massive stars can lose mass and angular momentum quickly 
enough to prevent the star reaching its critical rotation ve- 
locity on the main sequence. As a result, lower mass stars 
with lower mass loss rates would end their main sequence 
lives as more rapid rotators than higher mass stars. At 
the end of the main sequence, low mass rapidly rotating 
stars may reach the critical velocity as angular momen- 
tum is dredged up from the core. At the critical velocity, 
these stars would lose significant amounts of mass from the 
equator. This could explain why extremely rapidly rotating 
stars, such a s Be stars are not seen a mong the more mas- 
sive O stars (|Mevnet fc Maedeill2000l 'l. The interaction be- 
tween mass loss and angular momentum would also explain 
the relative frequency of Be stars in the LMC and SMC. 
At lower metallicity, the mass loss rates are typically lower, 
so stars in the LMC and SMC may be able to continue to 
rotate faster, and are he nce more likely to form Be stars 
ijMevnet fc Maedeij|200ol l 

A recent study by ICeorgy. Mevnet fc Maedeil (|201ll ') 
considers the effects of anisotropic winds on the evolution 
of a rotating star. Using the Roche approximation to deter- 
mine the shape of the surface equipotential, they find that 
as the rotation rate increases, more mass is lost from the 
pole. At high rotation rates, the difference in flux between 
the pole and equator becomes quite large, with the flux at 
the pole more than twice the flux at the equator when the 
rotation rate is greater than 0.8 Qc- Using a semi-analytical 
approach, when anisotropic effects are taken into account, 
the angular momentum loss changes by less than 4 % com- 
pared to isotropic winds for rotation rates up to 0.9 flc, and 
the subsequent effects on the stellar evolution are expected 
to be small. Using more realistic mass loss rates there is ba- 
sically no change in the angular momentum loss until Xc 
0.4. At this point, the star becomes sufficiently rapidly rotat- 
ing to significantly reduce the temperature of the equatorial 
regions, which pr oduces a strong equatorial mass l oss in the 
anisotropic case l|Georgv. Mevnet fc Maedeij|201l|) . 

An understanding of the latitudinal distribution of mass 
loss is also important for understanding the structure of the 
surrounding nebulae. The shape of the circumstellar neb- 
ulae may also have a n effect on supernovae light curves 
Ivan Marie et al.l HoiO). Early studies suggested that if the 
mass loss in a rotating star is primarily radial, the material 
carried away could lead to the formation of a wind com- 
pressed disk fWCD. lBiorkman fc CassincUi 1993). However, 
more recent research including the non-radial components of 
the wind have shown that this fiow inhibits the formation 
of a WCD, and may be m ore likely to produce a polar neb- 
ulae l|Owocki et al.(n"998l ). The effects of rotation on stellar 



lar density enhancements ( 
rial density enhancements ( 


Petrenz & Pulsl 2000f) or enuato- 


PeluDessv. Lamers & Vinklbood: 


iMadura. Owocki fc Feldmeieij 20071). Different mechanisms 



tion, as observations of B[e] stars ar e thought to be ex- 
plained by a dense equatorial w ind (jZickgraf et al.l Il985l : 
IZickgraf fc Schulte-Ladbec^ |l989!), while nebulae around 
LBVs are often bipolar, as is seen around rj Carinae. Asym- 
metric mass loss produced by rotation may be able to ex- 
plain the structure of some of these nebulae. 

Global effects of rotatio n on mass loss ha,ve been inves- 
tigated by several authors. iFriend fc AbbottI ([1986') found 
that the mass loss rate of a star at a given location in 
the HR diagram increases with increasing rotation rate . 
This effect has been confirmed observationally (IVardvall985l : 
INieuwenhuiizen fc do Jagcr 1988), and has also been found 
in models by Maedeii |1999). This relationship between mass 
loss and rotation rate can be described by a sim ple equation 
(eg, IBiorkman" fc Cassinellilll993l : lLangeJll998h : 



M{v) 



M(Vrot = 0) 



1 - V 



(1) 



for a star at constant temperature and luminosity. The ex- 
ponent ^ is usually taken to be 0.43. According to this 
equation, the enhancement of the mass loss rate due to ro- 
tation is modest, only a factor of 2 for Q = 0.8. Formally, 
this formula leads to infinite mass loss rates at critical rota- 
tion velocities, so this scaling law breaks down at very high 
rotation rates. It is this formula that is used to calculate 
the effects of rotation on mass loss in many stella r evolution 
codes (eg. lLangerlll998l : lMaeder fc Mevnetll2000l ) . 

The rotational enhancement of mass loss described 
above can be even more important in stars that are already 
close to the Eddington limit. In this case, for a star rotating 
at a given angular velocity Q, Equation [1] can be rewritten 
to include the effects of the Eddington luminosity as: 



M(n) 



(l-F)- 



M{Vrot = 0) 



n2 
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2-kGp„ 



may be important in different stages of stellar evolu- 



where pm is the average density inside a given surface, and 
a is one of the force multiplier parameters. In this case, 
the mass loss rate depends on both the rotation rate and 
the Eddington factor, a situation known as the fiT limit 
(|Maeder fc Mevnet|[200ol ). In extreme cases, with F > 0.639, 
even moderate rotation can cause the d enominator to van- 
ish, leading to infinite mass loss rates. iMaeder fc Mevne3 
(2000) show that the ratio M{Q.)/M{Q) can be quite mod- 
erate for stars up to 40 M©, but can become large for stars 
above 60 Mq , and can start to diverge for stars close to the 
Humphreys- Davidson limit. 

Despite recent progress, there is still much to be dis- 
covered on the interaction between stellar winds and rota- 
tion. Scaling laws, such as those used in Equations [1] and [2] 
are derived for stars at constant temperature and luminos- 
ity, but both of these quantities change as the rotation rate 
increases. Current ID and pseudo-2D models do not fully 
account for the distortion of the surface and variation in 
temperature, making realistic calculations of the anisotropy 
in the mass loss and angular momentum loss rates suspect. 
A more realistic calculation of these quantities could pro- 
duce significant discrepancies compared to current stellar 
evolution models. 

In this work, we calculate the stellar mass loss rates for 
fully 2D stellar structure calculations, addressing many of 
the issues discussed above. We consider two different mass 
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loss prescriptions, and compare the 2D results to ID calcula- 
tions and scaling laws. In this way, we can directly calculate 
the distribution of mass loss, and its effect on angular mo- 
mentum loss and the evolution of the star. 

In Section[2l we discuss the 2D stellar structure models 
used in this work. In Section [3] we compare our 2D calcu- 
lations to the mass loss rates and angular momentum loss 
rates in ID models. In Section |4] we discuss the effect of ro- 
tation on the global mass loss rates of stellar models, and 
present scaling relations derived at constant mass. Finally, 
we summarize our results and discuss future work in Section 

m 



2 MODELS 

2.1 Stellar structure 

The stellar models used here w ere calculated u sing t he 2D 
stellar structure code ROTORC jDeupred Il990l. Il995l). The 
code uses the OP AL opacities ( Iglesias fc Rogera 19961 ) and 
equation of state (|Rogers et all 19961 ). In this work, we con- 
sider ZAMS models of 20, 30 and 40 M©, with X = 0.7 and 
Z = 0.02. These models solve the conservation equations of 
mass, momentum, energy and hydrogen abundance, along 
with Poisson's equation for the gravitational potential on 
a two-dimensional finite difference grid with the fractional 
surface equatorial radius and the colatitude as the indepen- 
dent variables. This grid contains 581 radial zones and 10 
angular zones. Simultaneously with the conservation equa- 
tions, we solve the integral of the density over the model, 
and set the surface equatorial radius such that the result 
equals the total mass of the model. The surface is assumed 
to be an equipotential determined by the equatorial gravi- 
tational potential. The models, summarized in Table [1] are 
uniformly rotating ZAMS models, with velocities equal to 
0, 0.3, 0.5 and 0.7 n^. 



2.2 Including mass loss 

As ROTORC gives us the values of radius and effective temper- 
ature as a function of co-latitude, we can use these to calcu- 
late the local mass loss rate as a function of co-latitude in a 
rotating star. This has been done using two mass loss pre - 
scriptions: the numerical mass loss rates o f IVink et all ( |200lD . 
hereafter VOl, the origi nal theoretica l formu lation for radia- 
tively driven winds of ICastor et "al] (119751 ) (CAK) includ- 
ing th e finite disk correction, described by iKudritzki et all 
l|l989l ), hereafter K89. These two rates allow us to compare 
a theoretical formulation (K89) and the numerically derived 
rates commonly used in stellar evolution codes (VOl). 

In each angular zone, the local radius and effective tem- 
perature are calculated. These are used to calculate a local 
luminosity, which is defined as the luminosity of a spheri- 
cal star with the given radius and the effective temperature. 
Using these parameters, we can calculate a local Eddington 
factor, defined as 



(3) 



4-KcGM 

where ae is the electron scattering cross-section, which de- 
pends on the ionization fraction (and hence temperature 
and is hence allowed to vary as a function of co-latitude. 



If the vol mass loss rates are selected, the local Ed- 
dington factor is used to calculate the temperature of the 
bi-stability jump. The mass loss rate is then calculated in 
each angular zone using Equation 24 of VOl if the local 
effective temperature is above the temperature of the bi- 
stability jump, and their Equation 25 if the temperature is 
below the bi-stability jump. Below 12500 K, the VOl mass 
loss rates are undefined, and we neglect the mass loss from 
these regions. For the models considered here, all of the sur- 
face temperatures are greater than 25000 K, so this cutoff 
has no effect in our models. 

If the K89 mass loss rates are selected, the force mul- 
tiplier parameters (a, 5 and k) must also be determined 
locally. All of these parameters are temperature dependent, 
and we have u sed a parametr i zation based on the tabular 
values given in 'Leithere r et al (|2010l ). For k, we used data 
taken from their Tables 1 and 2 for solar metallicity stars 
of luminosity class V, and fit using a least squares analysis. 
This fit was implemented in our mass loss code, allowing k 
to be determined locally based on the effective temperature 
and luminosity. The values of a and 5 were set as a step 
function, with {a, 5) = (0.598,0.05) for Te// > 30000if and 
(0.607, 0.07) below this. For the K89 rates, we followed the 
method for calculating the mass loss rates outlined in K89, 
section 6, but we allow the luminosity, effective tempera- 
ture, and the force multipliers a, S and k to be functions of 
CO- latitude. 

The force multiplier parameter description of stel- 
lar winds is likely an over- simpli fi cation , and has been 
challenged by (for example) ILucvI l|2007h : iMiiller fc VinkI 
(2008) based on Monte Carlo simulations. A more detailed 
pa rametrization for the force multiplier has been described 
bv lKudritzkH (|2002l ). However, this new parameterization is 
much more complicated and must be solved iteratively. In 
the interest of computational efficiency, we have chosen to 
use the original force multiplier parameter description. 

Mass loss rates are typically given in terms of MQj/r~^, 
but this also has an implicit unit of per star. Since we are 
calculating local mass loss rates, we must scale the calcu- 
lated mass loss rates by the area of the zone. The calculated 
mass loss rates are multiplied by the time step to determine 
the total amount of mass lost, and this is then multiplied by 
the ratio of the area of the zone to the total area of a spher- 
ical star with the local radius to give the amount of mass 
lost from this zone. At each co-latitude, the mass is lost by 
removing surface zones until the required amount of mass 
has been lost. Finally, the total mass of the star is reduced 
by the amount of mass lost at this time step. The star is 
then allowed to relax back to an equilibrium configuration. 

As well as calculating mass loss in the zone-by-zone 
fashion described above, the code can also calculate a 
pseudo-ID mass loss rate, based on the luminosity and 
global effective temperature of the star. Calculating a 
pseudo-ID rate in this way allows us to make a more direct 
comparison with other ID calculations. We use the lumi- 
nosity and global effective temperature to calculate a single 
mass loss rate for the model. The total mass lost in this way 
is then divided among the angular zones weighted by area, 
the mass of the star is modified and the star is allowed to 
relax as described above. 

The evolution of a rotating star including mass loss pro- 
ceeds as follows. First, the star is allowed to evolve for some 
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Table 1. Model properties 



Mass 












AT 1 


L/Lq 


(M©) 


(km -1) 




(R©) 




(K) 


(K)" 




20 








5.835 


1.000 


34476 





42899 


20 


200 


0.3 


5.991 


0.969 


34090 


1161 


42313 


20 


375 


0.5 


6.437 


0.892 


33168 


3866 


41196 


20 


550 


0.7 


7.376 


0.770 


31802 


7899 


40122 


30 








7.310 


1.000 


39598 





117162 


30 


225 


0.3 


7.528 


0.967 


39195 


1470 


116564 


30 


400 


0.5 


8.038 


0.897 


38203 


4263 


113517 


40 








8.621 


1.000 


42959 





225732 


40 


250 


0.3 


8.893 


0.964 


42395 


1791 


222164 


40 


425 


0.5 


9.487 


0.897 


41336 


4711 


216740 



period of time. To lose mass, the evolution is halted and 
mass is lost as described above. Once the star has returned 
to an equilibrium configuration, the evolution is restarted, 
and continues until the next mass loss step. This process 
is repeated until the star reaches the desired evolutionary 
stage. In this paper, however, we will consider only zero-age 
main sequence (ZAMS) models. The effects of 2D mass loss 
on stellar evolution will be discussed in a future paper. 



3 RESULTS 

3.1 Mass Loss Rates 

The effects of rotation on mass loss rates can be seen in 
Figure [1] for the two considered mass loss rates. The VOl 
rates are always higher than the K89 rates, but the differ- 
ence decreases as the mass o f the mode l incre ases. This is 
consistent with the results of IVink et all (|2000l ). who found 
that their theoretical rates resolved the discrepancy between 
theoretical and observed mass loss rates for normal O stars. 
The observed rates are typically a factor of two higher than 
those predicted by radiation-driven wind theory, but the 
IVink et all (| 2000l) and vol rates are in good agreement with 
observations for these stars. 

In general, the mass loss rates decrease as the rotation 
rate increases. The exception is the K89 rate in the 30 Mq 
models, which increase slightly with increasing rotation rate. 
In this model, the effective temperature is close to 40 000 
K, where our parametrization of k has a discontinuity due 
to lines from CNO elements. Above 40 000 K, k increases 
sharply, and the resulting mass loss rates are higher than at 
temperatures below 40 000 K. In the 30 Mq models, the ef- 
fective temperatures are around 38 000-39 000 K, just below 
this gap. However, as the rotation rate increases, the pole 
becomes hotter, and temperatures in this region are above 
40 000 K. This local increase in the mass loss rate at the pole 
increases the overall mass loss rate for the star. The effect be- 
comes more pronounced as the rotation rate increases, and 
the mass loss rate increases with rotation. The 20 and 40 
Mq models are sufficiently far from this jump that the mass 
loss rates are not affected. The VOl models do not depend 
on our parameterization of the force multiplier parameters, 
and so this effect is not seen in these rates. 

Using our 2D models, we can also compare the 
anisotropy of the mass loss rates to that predicted using our 
pseudo-ID models. The difference between 2D and pseudo- 
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Figure 1. The mass loss rates as a function of rotation rate for 
the K89 rates (solid) and the VOl rates (dashed). From bottom 
to top, the rates are shown for 20, 30 and 40 Mq models. With 
the exception of the Kudritzki rates at 30 M©, the mass loss rates 
all decrease as the rotation rate increases. 



ID calculations is shown for the K89 rates in Figure [2] and 
for the vol rates in Figure [3| As expected, our 2D models 
lose more mass at the pole and less mass from the equator 
than the pseudo-ID models. The effects can be significant 
(nearly 0.1 dex at the pole) even for relatively slowly rotat- 
ing models (0.3 flc)- In Figure [2j the effect of the jump in k 
near 40 000 can be seen in the sharp decrease in mass loss 
rate near the equator of the 40 Mq model. Compare this 
with the 20 Mq model, where the temperature is always 
well below the jump, and the mass loss rate varies smoothly 
with co-latitude. 

Despite the differences in distribution of mass loss, the 
total integrated mass lost in the 2D and pseudo-ID calcula- 
tions is nearly identical. The total mass lost in each case is 
shown in Table [5| The difference in integrated mass lost is 
largest for models where the range of surface temperatures 
covers the jump in mass loss rates caused by k (30 and 40 
Mq models with the K89 rates) . The larger the temperature 
difference between the pole and the equator, the more zones 
are affected by the jump, and the bigger the difference in 
mass loss between the pseudo-ID and 2D models. With the 
vol rates, which are not affected by the jump in k, the dif- 
ference in total mass lost is typically less than 2 %, except 
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Figure 2. The variation with colatitude in the mass loss rates for 
the 20 M0 model (top) and the 40 M0 model (bottom). The log 
of the ratio of the 2D rates to the pseudo-ID rates are shown for 
the K89 mass loss rates at 0.3 f2c (solid) and at 0.7 Qc and 0.5 Qc 
(dashed) for the 20 and 40 Mq models respectively. Although the 
total mass lost is similar in the pseudo-ID and 2D calculations, 
more mass is lost from the pole in the 2D case, even at relatively 
slow rotation rates. 
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Figure 3. As for Figure (2] but using the VOl mass loss rates. 



in the most rapidly rotating models. In later phases of evo- 
lution, when the temperature of the star drops to the tem- 
perature of the bi-stability jump (22 500-26 COOK), a similar 
result can be expected in the VOl mass loss rates. 



Figure 4. Percent difference in angular momentum loss rate for 
the 2D and pseudo-ID mass loss rates in the 20 Mq model. Solid 
lines show the results using the K89 mass loss rates, while dashed 
lines show the results for the VOl rates. Shown are curves for 
models rotating at 550 kms~^ and 200 kms^^. 



predicted by the fully 2D mass loss rates, as shown in Ta- 
ble |3l The 2D model therefore loses less angular momentum 
through radiatively driven winds, and hence spin down more 
slowly than the ID models. As angular momentum is trans- 
ported from the core to the surface, a slower rate of angular 
momentum loss will result in either more rapidly rotating 
stars at the end of the main sequence, or more mass lost at 
the equator as the star reaches its critical velocity. 

As expected, as the rotation rate of a model increases, 
the rate of angular momentum loss increases as well, as il- 
lustrated in Figure [5] However, at higher rotation rates the 
total angular momentum of the models also increases. For 
each model and mass loss rate, we calculated the ratio of 
the local rate of angular momentum loss to the total angular 
momentum loss {L{6)/L). For a given model, this quantity 
is relatively constant as rotation rate increases, as shown 
in Figure |6] for the 20 Mq model. This is true for any given 
model, but different masses and mass loss prescriptions scale 
differently, making it difficult to derive a simple scaling law 
to describe the angular momentum loss. In more rapidly ro- 
tating and more massive models, there is more variation in 
mass loss near the equator, so the fractional rate of angular 
momentum loss does decrease as the rotation rate increases. 
This is particularly true for the VOl models, while the K89 
models show a more consistent scaling. 



3.2 Angular Momentum Loss 

Our technique for modelling the mass loss in 2D also allows 
us to directly determine not only the total rate of angu- 
lar momentum loss, but also the rate of angular momen- 
tum loss from each surface zone. This can be done for both 
the pseudo-ID and the 2D mass loss rates described in sec- 
tion 12.21 The surface variation in angular momentum loss 
is shown in Figure [4] for the 20 Mq models rotating at 200 
and 550 km s~^. The overall rate of angular momentum loss 
predicted by the pseudo-ID mass loss is higher than that 



4 DISCUSSION 
4.1 Global mass loss 

The current paradigm for stellar mass loss uses Equation [1] 
to scale mass loss rates with stellar rotation for stars at con- 
stant luminosity and effective temperature. As discussed in 
Section [Tl this effect has been seen obs ervationally ( Vardyal 
Il985l : iNieuwenhuiizen Sz de Jage"illl988l ). although the mag- 
nitude of the effect appears to be somewhat uncertain, rang- 
ing from a few tens of percent to two orders of magnitude. 
However, stellar models use mass as the independent vari- 
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Table 2. Total mass (Mq) lost over 1000 years 









K89 






vol 




Mass 


Veq 


pseudo-ID 


2D 


% difference 


pseudo-ID 


2D 


% difference 


20 


200 


1.0463e-6 


1.0436e-6 


0.26 


2.4137e-5 


2.4148e-5 


0.045 


20 


375 


9.7383e-7 


9.8561e-7 


1.2 


2.1206e-5 


2.1598e-5 


1.8 


20 


550 


9.1649e-7 


9.5660e-7 


4.2 


1.7708e-5 


1.9337e-5 


8.4 


30 


225 


1.9915e-5 


2.3663e-5 


16 


1.6802e-4 


1.6802e-4 


0.0012 


30 


400 


1.8372e-5 


2.7743e-5 


34 


1.5354e-4 


1.5473e-4 


0.77 


40 


250 


2.6088e-4 


2.5734e-4 


1.4 


5.0327e-4 


5.0317e-4 


0.020 


40 


425 


2.4987e-4 


2.1379e-4 


17 


4.7070e-4 


4.7230e-4 


0.34 



Table 3. Total angular momentum lost over 1000 years 









K89 






vol 




Mass 


Veq 


pscudo-lD 


2D 


% difference 


pseudo-ID 


2D 


% difference 


20 


200 


3.06096c47 


2.93396e47 


4.3 


6.35610e48 


6.19461e48 


2.6 


20 


375 


5.54387e47 


4.90389c47 


13 


1.10423e49 
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Figure 5. Rate of angular momentum loss calculated in 2D for 
a 20 Mq model using the K89 (top) and VOl (bottom) mass loss 
rates. The three curves show the variation in the rate of angular 
momentum loss across the surface for models rotating at 0.3 Qc 
(solid), 0.5 (dashed), and 0.7 Qc (dot-dashed). 



Figure 6. Local rate of angular momentum loss years scaled 
by the total angular momentum loss of the model for 20 M© 
models rotating at 0.3 flc (solid), 0.5 He (dashed), and 0.7 Qc 
(dot-dashed). Angular momentum loss was calculated using the 
K89 (top) and VOl (bottom) mass loss rates. 



able, and the luminosity and effective temperature are al- 
lowed to vary. 

Calculations of fully 2D rotating stars indicate that the 
situation may be more complicated, and Equation [T] may be 
used incorrectly in some situations. We find that for models 
at constant mass, the overall mass loss rates decrease as the 
rotation rate increases, as shown in Figure[7]for 20 Mq mod- 
els. As the rotation rate increases, the effective temperature 
and luminosity also decrease, as shown in Table [1] which 
leads to lower mass loss rates in all cases. The decrease is 
about 0.06 dex for the K89 rates, and about 0.15 dex for the 
vol rates. At 30 and 40 Mq, the decrease is between 0.03 



- 0.05 dex for all three rates, although these models do not 
extend to as rapid rotation as the 20Mq models. 

As discussed above, the scaling law in Eqn [1] applies to 
stars at constant luminosity and effective temperature, but 
is often applied to models with the same mass at different 
rotation rates even though the rotation changes Teff and L. 
The mass loss rates discussed here are held at constant mass, 
while the effective temperature and luminosity are allowed 
to vary. This makes more sense from a stellar modelling 
stand point as in most evolution codes, models are made at 
a given mass, and the effective temperature and luminosity 
are allowed to vary as the rotation rate changes. For those 
who wish to implement a mass dependent scaling law, we 
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Figure 7. The 2D mass loss rates for a 20 Mq model as a function 
of rotation rate (solid curves) for the K89 (bottom) and VOl 
(top) mass loss prescriptions. The mass loss rates predicted by 
the scaling law given in Eqn[T]are shown by the dashed curves. 



have found linear fits to our 2D calculations. For the VOl 
models, the mass loss rate can be fit linearly in the mass 
and the square of the ratio of the surface equatorial rotation 
velocity to the critical velocity {{v/vcrit)^)- 

M{v) = M{v = 0) ~0.464*{v/vcr^tf +0.007157{M/Mq).{4) 

The K89 rates are more complicated. At any given mass, the 
mass loss rate can be fit linearly as a function of (v/vcrit)^, 
but the slope of the line varies with mass. With only three 
points, it is difficult to produce a reliable fit, although the 
result is clearly not linear. We present here fits at each mass 
individually: 



M{v) = M{v = 0) + a{v/vcrit) 



(5) 



where a = -0.148 at 20 Mq, 
40 M0. 



-0.221 at 30 Mq and -0.123 at 



4.2 Anisotropic mass loss 



As discussed in Section [STTJ one of the consequences of rota- 
tion is to produce a variation in temperature across the sur- 
face of the star, which causes the mass loss rates to vary as 
well. Based on the CAK mass loss formalism, the mass loss 
rate at any given colatitude in a rotating star can be written 
relative to the polar values. Mo, Fo and go = GM/B?: 



(6) 



where a is taken to be independent of latitude and g^ff is 
the centrifugally reduced surface gravity 



M{e) 


\m] 


l/a 


'geffiey 




Mo 


Fo 




9° 





(7) 



where D.'^ = Vd^i^/GM fOwocki et all (Il996l . ll998t ). Combin- 
ing this equation with the von ZeipeT(|l924 ') gravity darken- 
ing shows that M{6) oc ffe ffjd) , and hence th e mass loss rate 
decreases towards the equator lOwocki et al.l l| 19961 ). This is 
expected to produce bipolar outfiows around rotating stars. 



Figure 8. The relative difference between the 2D calculation and 
the ID prediction for mass loss rate from a 20 Mq model rotat- 
ing at 200 (top) and 550 kms~^ (bottom panel)as a function of 
colatitude for the K89 (solid) and VOl (dashed) mass loss rates. 
The ID prediction is scaled such that the total mass loss is the 
same. 



Using our 2D stellar models, we can test the efficacy of this 
scaling law approximation. 

The difference between the ID predictions and the 2D 
calculations for mass loss rate across the surface of a rotating 
star is shown in Figure [8] for a 20 Mq model rotating at 200 
and 550 fems~^and in Figure [9] for a 40 Mq model. We 
have scaled the mass loss rates predicted by Equation |B] so 
that the total mass loss is the same. In all cases, the ID 
scaling law over predicts the amount of mass loss from the 
equatorial regions, and under predicts the mass loss from 
the polar regions. The differences between the VOl models 
and the ID predictions are typically larger than the K89 
rates by a factor of 2. In the most rapidly rotating model 
(20 Mq, V — 550 kms~^), the VOl mass loss rates are 100 
% smaller than the ID predictions at the equator. For most 
models, typical variations for all mass loss rates are between 
10-20 % larger at the pole and smaller at the equator. 

Most of the difference between the ID and 2D models 
can be accounted for by differences between the tem pera- 
tures predicted by von Zeipel's law Ivon Zeipel (Il924l ) and 
the temperatures calculate d in our models, wh i ch are typi- 
cally in the range of 10 % iLovekin fc Deupred l|2006l ). The 
remaining difference arises from the assumptions used in de- 
riving Equation!?] In reality, the force multiplier parameters, 
a, 5 and k are all functions of temperature and luminosity. 
In Equation [7] however, a is taken to be independent of lat- 
itude, and only the effective gravity varies across the sur- 
face. In our models, a is nearly independent of temperature, 
and so we also assume a to be constant. The second pa- 
rameter, 5, is also approximately constant over the range of 
temperatures considered here. However, the third line force 
parameter, k, is strongly dependent on temperature. This 
dependence on temperature, which is not accounted for in 
the ID scaling relation, will increase the variation between 
pole and equator in fully 2D mass loss rates. 
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Figure 9. The relative difference between tlie 2D calculation and 
the ID prediction for mass loss rate from a 40 Mq model rotating 
at 250 (top panel) and 400 ferns" ^ (bottom panel) as a function 
of colatitude for the K89 (solid) and VOl (dashed) mass loss rates. 
The ID prediction is scaled such that the total mass loss is the 
same in both cases. 



5 CONCLUSIONS 

Two difTerent prescriptions for radiatively driven mass loss 
rates we r e add e d to the 2D stellar evolution code RDTDRC 
iDeupred (|l99(]l . 1 19951 ). We used (1) a modification of the 
original theoretical d escription which inclu des the finite 
disk approximation bv lKudritzki et all ()l989l ) an d (2) mass 
loss ra tes based on Monte Carlo simulations bv IVink et all 
l|200ll ). Either of these mass loss prescriptions can be used 
to calculate the mass loss using the local properties across 
the surface of a rotating star, or using the global properties 
to give a pseudo-ID calculation of the mass loss. 

We considered the total mass lost from rotating stars 
using our models. We found that as the rotation rate in- 
creases, the mass loss rate decreases for any given mass. 
As rotation increases the effective temperature and lumi- 
nosity of the star decreases, which decreases the mass loss 
rate. This effect goes in the opposite direction of the com- 
monly used scaling relation given in Equation [T] However, 
this equation is derived, both theoretically and observation- 
ally, for stars held at constant temperature and luminosity, 
which is not the case when comparing stars of the same 
mass and increasing rotation rate. In Equations U and [S] we 
provide scaling relations for mass loss rates as a function of 
rotation rate for stars at constant mass. 

Using our 2D models, we investigated the variation in 
mass loss rate across the surface of a rotating star. A di- 
rect comparison of our pseudo-ID and 2D calculations show 
that the total angle integrated amount of mass lost is not 
significantly affected, but a full 2D calculation dramatically 
changes the distribution of mass lost. As our models do not 
assume von Zeipel's law holds, wc can provide an indepen- 
dent verification of the scaling relation b as ed o n the CAK 
theory and von Zeipel's law t\von Zeipeli il924i ) commonly 
used in ID models. We have found that the von Zeipel mod- 
els predict a smaller range of mass loss rates from the pole 
to equator than our 2D models, with typical differences of 
10-20 % at the pole and equator. The difference between 



the 2D and pseudo-ID model can mostly be explained by 
differences between our surface temperature structure and 
that predicted by von Zeipel's law, with the remaining dif- 
ferences a result of allowing the force multiplier parameters 
to vary with co-latitude in our 2D models. 

The change in distribution of the mass loss also affects 
the angular momentum loss from these models. In general, 
the 2D models predict a lower rate of mass loss at the equa- 
tor than the pseudo-ID models, for which the mass loss rate 
is constant across the surface. Although the total rate of 
angular momentum loss increases as the rotation rate in- 
creases, the fraction of the total angular momentum lost is 
approximately constant for a given mass and mass loss pre- 
scription. However, the scaling changes dramatically with 
mass and mass loss prescription, so it does not seem that a 
straightforward scaling relation can be derived. 

As the 2D models presented here lose less angular mo- 
mentum, they should be expected to spin down more slowly 
than ID models. This could produce one of two effects. The 
first possibility is that these models will lose more mass at 
the pole and remain rapidly rotating through later stages of 
evolution. The second possibility is that these models would 
maintain a rotation rate high enough to reach the critical ve- 
locity sooner in the evolution and lose a significant amount 
of mass (and hence angular momentum) at the equator. In 
a future paper, we will investigate how the 2D mass and 
angular momentum loss affect the evolution of these stars. 
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